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ABSTRACT 


Twenty-four  high-resolution  FMCW  soundings  taken  over  an  east-west 
1900-km  forward-propagation  path  show  unexpectedly  complex  mode  struc¬ 
ture  near  sunset.  Because  of  the  presence  of  1-  and  2-hop  sporadic-E 
modes,  it  was  suggested  that  multiple  reflections  from  sporadic-E  and 
F-layers — such  as  "m"  and  "n"  modes--ac.counted  for  much  of  the  unexplained 
structure . 

This  work  represents  an  attempt  to  see  how  well  these  "mixed"  modes 
can  be  understood  by  employing  computer  ray-tracing  procedures  developed 
by  Dr.  T.  A.  Croft  at  Stanford  Electronics  Laboratories. 

A  simple  p-Chapman  layer,  with  its  parameters  adjusted  to  give 
synthetic  ionograms  that  matched  the  experimental  1-  and  2-hop  F-layer 
modes,  proved  to  be  an  adequate  model  for  the  F-layer  profile  over  the 
path.  Reference  to  the  calibrated  experimental  ionograms  yielded  values 
for  the  sporadic-E  reflection  heights  for  each  time  of  day.  Assuming  no 
horizontal  tilts  in  layer  density  and  ignoring  the  earth's  magnetic  field, 
synthetic  ionograms  involving  multiple  F-layer  and  sporadic-E  reflections 
were  then  computed . 

The  model-experiment  comparisons  show  that  surprisingly  accurate 
estimations  of  both  1-  and  2-hop  and  mixed  modes  are  possible,  using  only 
relatively  simple  computational  procedures.  Although  similar  work  has 
been  done  in  this  area,  the  present  study  is  novel  in  its  use  of  calibra¬ 
ted  oblique  ionograms  for  model  derivation.  The  work  should  lead  to  a 
better  understanding  of  both  HF  forward  and  backscatter  propagation 
studies  which  are  plagued  with  multiple,  interlayer  propagation  modes. 

A  further  extension  of  the  modeling  technique  would  include  the 
insertion  of  additional  electrons  in  the  ionospheric  profile — E-,  F  -, 
and  F-layers — as  becomes  necessary.  Horizontal  tilts  in  electron 

a 

density  could  also  be  modeled  by  use  of  a  slow1 v-vary ing  profile  or  by 
insertion  of  two  or  more  different  discrete  models  along  the  propagation 
path . 
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s 
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Kf 

N  (h) 


constant  used  in  computing  a  discrete  tilt 
constant  used  in  computing  a  discrete  tilt 
constant  used  in  computing  a  discrete  tilt 
equivalent  E^-layer  ground  distance 
equivalent  F-layer  ground  distance 
a  computer  solution  for  ground  distance 
path  length,  1881.27  km 
frequency 

profile  shape  factor 

discrete  F-layer  virtual  height  at  point  "a" 
discrete  F-layer  virtual  height  at  point  "b" 
discrete  F-layer  virtual  height  at  point  "c" 
virtual  height  of  sporadic-E  layer 
virtual  height  of  F-layer 
height  of  maximum  ionization 
height 

scale  height  that  appears  in  p-Chapman  equation  for 
density 

scaling  factor  of  model's  density  function 
constant  used  in  computing  a  variable  tilt 
constant  used  in  computing  a  variable  tilt 
frequency  scaling  factor,  A 
electron  density  as  function  of  height 


N 


N' 
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maximum  electron  density 

2 

revised  maximum  density.  =  K  N 

J  ’  f  o 

specified  constant  used  for  calculating 
specified  constant  used  for  calculating 
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specified  constant  used  for  calculating  D^, 
specified  constant  used  for  calculating  TQ 
specified  constant  used  for  calculating 
mean  earth  radius,  6370  km 
group  time  delay 

equivalent  group  time  delay  for  a  single  hop  to  the 
F-layer 

equivalent  group  time  delay  for  a  single  hop  to  the 

E  -layer 
s 

(height  -  h  )/H 
ms 

take-off  (elevation)  angle 

(actual  range  -  starting  range )/1881 .27,  variable  used 
in  modeling  a  variable  tilt 

index  of  refraction  of  ionosphere 


solar  zenith  angle 
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I  INTRODUCTION 


A .  Purpose 

The  purpose  of  the  work  reported  here  was  to  determine  the  origin 
of  a  set  of  unexplained  modes  appearing  on  a  series  of  oblique  sounding 
records  of  high-frequency  (HF)  radio  signals  taken  over  an  E-W  1900-km 
path  between  Stanford  University,  California,  and  Lubbock,  Texas.  It 
was  suggested  that  the  mysterious  modes  were  the  result  of  interlayer 
propagation  between  sporadic-E  (E^)  and  F-layers.  The  author  attempted 
to  verify  this  suggestion  by  employing  currently  available  computer  ray¬ 
tracing  techniques  to  aid  in  the  modeling  of  several  multiple-layer 
propagation  modes,  which  were  then  compared  with  the  experimental  records. 

B .  Background 

In  the  late  1950s,  researchers  began  to  consider  seriously  the 
different  possible  modes  and  corresponding  angles  of  arrival  that  might 
occur  when  HF  radio  signals  are  transmitted  over  a  given  path. 

In  1953,  Wilkins  and  Kift  [Ref  .  1],  working  at  the  Department  of 
Scientific  and  Industrial  Research  (DSIR)  Radio  Research  Station  in 
Slough,  England,  measured  angles  of  elevation  of  incoming  radiation  from 
pulse  and  telegraph  signals  originating  in  Negombo,  Ceylon,  and  Kirkee, 
India,  and  from  telegraph  signals  originating  in  New  Delhi.  For  frequen¬ 
cies  between  15  and  19  MHz,  energy  was  received  most  often  and  most 
strongly  in  the  elevation  range  7  ±  2  degrees.  In  the  summer  of  1954, 
signals  from  Kirkee  showed  a  decrease  in  elevation  at  the  receiver  to 
5  ±  2  degrees;  this  circumstance,  correlated  with  backscatter  observations 

at  both  Ceylon  and  Slough,  suggests  that  th~  summer  increase  in  E  -layer 

s 

ionization  causes  the  decrease  in  received  angles  of  elevation. 

In  October  1956,  Silberstein  [Ref.  2]  at  the  National  Bureau  of 
Standards  (NBS)  in  Boulder,  Colorado,  ran  a  long-distance  pulse-propaga¬ 
tion  experiment  at  20.1  MHz  between  Sterling,  Virginia,  and  Maui,  Hawaii. 
Simultaneous  oblique  ionograms  were  made  along  the  path  from  Sterling  to 
Boulder,  and  vertical  incidence  records  were  also  obtained  at  Sterling. 
Silberstein 's  results  showed  mode  structures  that  differed  greatly  from 
one  day  to  the  next,  indicating  that  a  long  path  is  very  sensitive  to 
ionospheric  conditions. 
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In  1959,  Kift  [Ref.  3]  remarked  that  .  .as  early  as  1938, 
Millington  .  .  .  pointed  out  that  long-distance  paths  would  consist  of 
hops  of  unequal  length,  and  that  a  knowledge  of  the  distribution  of 
ionization  at  every  point  along  the  route  was  desirable  in  solving  the 
problem  of  how  energy  travels  from  one  point  on  the  earth's  surface  to 
another."  Kift  assumed  that,  in  all  thick  layers,  the  distribution  of 
electron  density  vs  height  followed  a  parabolic  function,  and  that  no 
lateral  tilts  of  the  layers  were  present.  It  was  assumed  that  reflection 
from  the  E^-layer  was  analogous  to  that  from  a  mirror.  Kift's  treatment 
of  the  problem  is  based  on  the  parabolic-layer  transmission  equations  of 
Appleton  and  Beynon  [Ref.  4],  and  uses  ray  theory  throughout.  He 
examined  two  long  paths  and  found  that,  when  the  effect  of  tropical  E 

s 

ionization  was  included,  he  could  predict  time  delays  and  angles  of 
arrival  for  the  various  modes  with  some  confidence.  He  used  a  fixed- 
frequency  pulse  sounder  system,  which  yielded  measurements  of  angles  of 
arrival  and  propagation  times  for  the  various  modes.  He  then  compared 
the  results  with  those  predicted  on  the  basis  of  a  model  assuming  a 
parabolic  layer. 

Investigators  at  Stanford  Research  Institute  (SRI),  under  the 
supervision  of  R.  Vincent,  did  some  extensive  research  titled  "HF  Commu¬ 
nication  Effects,"  which  was  finished  by  1965;  Ref.  5  gives  a  summary  of 
their  activity.  Several  oblique  HF  communication  paths,  including 
midlatitude  and  transauroral,  were  studied.  A  primary  interest  was 
given  to  what  was  termed  the  "Channel  Scattering  Function, "  with  which 
the  HF  channel  is  characterized  by  a  3-dimensional  plot  of  amplitude  vs 
time  delay  and  frequency.  Individual  modes  (both  simple  and  complex) 
are  separated  by  their  time  of  travel  over  the  path,  as  well  as  by  their 
characteristic  doppler  shift  due  to  layer  motion  and  changing  density. 

The  phase-stable  system  used  for  these  measurements  is  discussed  in 
Ref.  6.  The  application  of  the  phase-stable  system  to  the  measurement 
of  frequency  dispersion  and  doppler  shift  over  both  a  transauroral  and  a 
midlatitude  path  is  outlined  in  Refs.  7  and  8.  Simultaneous  phase-stable 
transmissions  on  7.366  MHz,  paralleled  by  step-frequency  oblique  incidence 
sounders,  enabled  the  SRI  workers  to  identify  complex  modes  involving 
interlayer  propagation  between  sporadic-E  and  F-layers — such  as  "M"  and 
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to  Plot  f  E  „„  investigators  used  their  sounder  data 
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own  that  illustrates  the  relative  placement  of  the  various 

modes  on  an  exampie  oblique  ionogram.  (Several  plots  are  then  presented 
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,  nlotted  vs  time  of  day,  for 

that  show  "accumulated  percentage  by  mode  plott 

several  different  paths  and  seasons.)  It  was  found,  for  examp  e, 

rincipal  daytime  modes,  over  a  path  item  Palo  Alto,  California  to 
Zl  Greenland,  ion  winter  conditions,  .ere  in  decreasing  order  - 

significance :  the  2F,  H,  ».  and  IF.  the 

Quite  recently  (1967),  P.  L.  George  [Ref.  13],  «o 

Australian  Department  of  Supply,  conducted  a  series  of  HP 

experiments  over  a  1900-km  path  in  Australia,  located  near  ? 

Magnetic  Meridian.  He  ohtained  Doth  obligee  soundtngs  -* 

soundings  at  the  midpoint  of  the  path.  He  assumed  a  par  hoi  c 

model,  which  .as  determined  from  the  midpoint  records  taken 

•  S  His  model  contained  a  sporadic-E-layer,  an  E-layer,  an 
frequencies.  His  moae  E  ) 

F  -layer  He  then  synthesized  an  ionogram  containing  , 

2  I  IE  2F  3F  4F,  and  E  -F  modes.  The  synthesized  ionograms 

2Es’  S’  ud  ’  tch  to  the  corresponding  experimental  records  for  IE,,  ^ 

r:r.r: ::  :i „ «. „  reproduced 

George ’  s  report.  There  also  seems  little  douht  about  the  identification 


**  FREQUENCY  (MHZ) 

SAMP1E  OF  a  —FEED  OBEIQDB  iOHOGRAM  ^P—  0H»  EXPERiMEHTAE 
(ECORI)  TAKEN  OVER  A  1900-km  PATH, 


of  the  3F  and  4F  modes,  although  there  was  a  discrepancy  of  a  few  tenths 
of  a  millisecond  between  the  experimental  results  and  the  predictions 
based  on  the  model.  It  is  difficult  to  predict  the  location  of  George's 
VF  mo<^e>  because  the  record  was  highly  cluttered  in  the  region  between 
the  IF  and  2F  traces . 

In  summary,  it  appears  that  Dayharsh  and  George  have  made  the 
greatest  contributions  toward  explaining  complex  mode  structure  by  use 
of  computer  synthesis  procedures.  Although  the  paths  that  Dayharsh 
studied  were  longer  and  more  disturbed  than  was  that  studied  by  George, 
Dayharsh  had  access  to  more  vertical  soundings  from  which  he  could  model 
a  tilted  ionosphere.  Therefore,  the  reliability  of  the  two  techniques 
is  probably  comparable.  The  major  problem  in  any  analysis  of  this  sort 
is  that  the  ionospheric  profile  actually  remains  unknown,  and  a  good 
guess  is  made  on  the  basis  of  the  behavior  of  HF  propagation  through  the 
various  layers  of  electron  density .  The  density  profile  is  therefore 
derived  by  an  indirect  method,  over  a  specific  point  along,  or  near,  the 
path,  and  is  then  assumed  to  vary  smoothly  (or  perhaps  stepwise)  along 
the  path. 

E.  M.  Young,  one  of  the  SRI  workers  involved  in  the  study  discussed 
above,  notes  [Ref.  14,  p.  16],  ".  .  .a  theory  that  would  permit  predic¬ 

tion  of  complex  modes  in  satisfactory  agreement  with  measured  data  has 
not  yet  been  developed."  Within  the  above  context,  the  present  work  will 
provide  still  another  contribution  in  the  area  of  complex-mode  synthesis 
by  computer  ray-tracing  techniques.  However,  the  primary  differences  in 
the  present  work  are  (a)  no  vertical  soundings  were  available  with  which 
one  could  determine  ionospheric  structure  along  the  path;  therefore  (b) 
a  p -Chapman  ionosphere  is  assumed  for  the  F-layer.  Thus,  only  a  sequence 
of  calibrated  oblique  ionograms  is  used  in  the  present  work  to  perform 
the  complex— mode  analysis;  the  ionospheric  models  are  derived  by  making 
an  educated  guess  at  a  first  model,  and  then  modifying  parameters  until 
a  suitable  final  model  is  obtained.  Tilts  could  be  modeled  only  by  exam¬ 
ining  model  variation  with  time  of  day,  because  each  model  is  obtained 
through  matches  to  a  sequence  of  IF  oblique  modes. 
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Approach  Used  in  the  Present  Work 


The  aforementioned  oblique  soundings  taken  over  an  E-W  1900-km  path 
between  Stanford  University,  California,  and  Lubbock,  Texas,  showed  some 
complicated  mode  structure  toward  the  latter  part  of  the  day .  At  the 
same  time,  a  sporadic-E  layer  became  evident,  and  it  was  conjectured 
that  the  "in-between"  modes  were  those  resulting  from  multiple,  inter¬ 
layer  propagation  among  the  E  -  and  F-layers.  Midpoint  vertical 
soundings  were  not  available;  however,  the  oblique  ionograms  were 
calibrated  in  time  delay. 

Fortunately,  computer  ray-tracing  techniques  had  already  been 
developed  at  Stanford  by  Dr.  T.  A.  Croft.  These  techniques  had  been 
applied  to  the  tracing  of  rays  through  model  ionospheres.  From  these 
ray  tracings,  oblique  ionograms  could  be  constructed  [Ref.  15].  The 
techniques  had  not  been  adapted  for  use  in  the  construction  of  models 
of  multiple  layer  mode  ray  sets;  however,  after  a  short  inquiry,  it  was 
deemed  feasible  to  attempt  such  a  construction. 

It  was  decided  to  apply  ray-tracing  techniques,  assuming  a  p-Chap- 
man  layer  [Ref.  16,  p.  17]  as  a  model  for  the  F-layer,  in  order  to: 

(1)  see  whether  by  this  procedure,  a  sufficiently  accurate  model 
could  be  obtained  for  the  F-layer  without  having  to  include 
the  effect  of  horizontal  tilts  in  the  electron  density  profile; 
and 

(2)  determine  whether  the  various  multi-hop  modes  could  be  sucess- 
fully  synthesized.  If  so,  the  result  could  explain  the 
mysterious  structure  that  gives  rise  to  the  unexplained  modes 
appearing  over  the  transmission  path. 

This  report  follows  through  the  process  by  which  the  synthesized 
ionograms  were  obtained  and  compared  with  the  experimentally  obtained 
ionograms.  The  results  are  also  compared  with  those  derived  on  the 
basis  of  a  simpler,  mirror-analogy  ionospheric  model.  The  sources  and 
magnitudes  of  errors  are  then  estimated.  The  effect  of  tilts  in  electron 
density  profiles  (and  how  to  model  them)  is  discussed  with  regard  to 
possible  future  work  in  mode  prediction  studies. 
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II  THE  EXPERIMENT 


A  series  of  PMCW  oblique  ionogram  records--cal ibrated  both  in  time 
delay  and  in  frequency — was  taken  over  a  path  between  Stanford,  California, 
and  Lubbock,  Texas,  on  Feb.  27,  1967;  the  path  length  was  1881.27  km; 
samples  of  these  records  are  presented  in  Figs.  1-24,  and  will  be 
discussed  subsequently.  Much  greater  clarity  and  larger  signal-to-noise 
ratio  were  achieved  on  these  records  than  was  ever  before  possible — 
especially  at  the  lower  frequencies — using  relatively  low  power.  This 
improvement  over  previous  ionogram  records  is  attributable  to  the  FMCW 
technique  developed  by  Drs .  R.  B.  Fenwick  and  G.  H.  Barry  while  working 
with  some  of  the  staff  at  the  Radioscience  Laboratory  at  Stanford  Univer¬ 
sity  . 

The  equipment  used  has  been  described  elsewhere  [Refs.  17-20],  The 
time-delay  calibration  was  achieved  by  carefully  calibrating  two  Hewlett- 
Packard  (HP)  cesium-beam  standards,  one  of  which  was  subsequently  flown 
to  the  field  site  in  Lubbock.  These  standards  yielded  a  time  calibration 
of  the  system  starting  pulse  and  constituted  the  frequency  standard  from 
which  each  sweep  system  derived  its  calibrated  input  signal. 

One  hundred  records  were  taken  every  hour,  over  a  24-hour  period. 

A  motion  picture  was  made  to  show  a  running  sequence  of  the  records 
throughout  the  day.  Time  and  frequency  scales  and  notices  of  time  of 
day  were  also  included.  This  picture  was  run  through  a  number  of  times; 
several  interesting  features  in  the  ionograms  that  changed  with  time  of 
day  became  apparent. 

One-  and  two-hop  F-layer  returns  were  noted  throughout  the  day  into 
early  nighttime.  Vertical  time-delay  extensions  on  the  IF  traces  noted 
later  in  the  evening  can  probably  be  attributed  to  multiple  reflections 
from  ionospheric  irregularities  (spread-F)  [Ref.  7,  p.  153].  Evidence 
of  ground  forward  scatter — reflections  from  irregularities  on  the  ground 
located  at  arbitrary  points  along  the  path — causing  a  severe  vertical 
broadening  of  the  2F  and  3F  traces,  was  present  throughout  the  day. 

Some  interesting  traveling  ionospheric  disturbances  (TIDs)  appeared 
several  times  in  the  record  sequences;  these  have  subsequently  been 
investigated  by  Dr.  T.  A.  Croft  at  Stanford. 
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The  most  notable  feature  of  the  records  was  the  appearance  of 
several  "in-between"  modes — modes  appearing  between  IF,  2F,  and  3F 
traces — beginning  in  the  later  part  of  the  day .  As  these  modes  became 
more  pronounced,  it  was  also  noted  that  a  sporadic-E-layer  had  formed, 
yielding  1-  and  2-hop  Eg  modes — which  may  be  identified  from  the  time- 
delay  vs  frequency  data  as  appearing  as  nearly  straight  lines  with  very 
high  maximum  usable  frequencies  (MUFs).  The  fact  that  the  lE^  trace  did 
not  extend  into  the  lower  frequencies  earlier  in  the  day  suggests  high 
D-layer  absorption,  coupled  with  low-gain,  low-angle  antenna  patterns  at 
low  frequencies.  Later  in  the  evening,  the  1E^  traces  became  much  more 
pronounced,  as  one  would  expect . 

To  investigate  the  causes  of  these  in-between  modes,  24  records  at 
various  times  throughout  the  day,  with  the  majority  at  sunset,  were 
selected  (Figs.  1-24);  these  photographs  were  enlarged  to  approximately 
8  by  10  inches,  a  convenient  size  for  handling  and  scaling.  Calibrated 
scale  underlays  on  which  the  photos  could  be  laid  over  a  light  box  were 
constructed.  The  calibration  was  achieved  by  noting  the  time-delay  tick 
marks  on  the  photos  and  by  deducing  the  frequency  scale  from  the  horizon¬ 
tal  end-points  of  the  records  that  contained  the  linear  frequency  sweep 
with  known  start  and  stop  frequencies.  The  underlays  were  used  to  plot 
the  various  models  derived  as  explained  below. 

Since  the  unexplained  modes  were  not  present  when  the  E^  reflections 
were  also  absent — as  Figs.  1-24  show  when  early  and  late  times  of  the 
day  are  compared — it  was  suggested  that  the  in-between  modes  were  a 
result  of  multiple,  alternate  reflections  from  the  E  -  and  F-layers  (as 
previously  mentioned).  There  seems  to  be  no  indication  of  the  presence 
of  a  normal  E-layer,  since  the  characteristic  S-turn  transition  between 
the  lower  IF  and  upper  IE  traces  was  never  evident.  The  S-turn  maj  have 
occurred  at  some  frequency  lower  than  4  MHz;  it  therefore  would  not  be 
visible  within  the  data  taken.  However,  if  this  did  occur,  it  is  not 
likely  that  a  layer  of  such  low  density  would  seriously  perturb  the  mode 
structure,  and  it  could  therefore  be  ignored. 
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Fig.  1.  EXPERIMENTAL  CALIBRATED  OBLIQUE  IONOGRAM 
MLT  =  Midpath  local  time. 
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Fig.  8.  EXPERIMENTAL  CALIBRATED  OBLIQUE  IONOGRAM 
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EXPERIMENTAL  CALIBRATED  OBLIQUE  IONOGRAM  11 


EXPERIMENTAL  CALIBRATED  OBLIQUE  IONOGRAM  12 


EXPERIMENTAL  CALIBRATED  OBLIQUE  I ON OG RAM  13 


EXPERIMENTAL  CALIBRATED  OBLIQUE  IONOGRAM  16 
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EXPERIMENTAL  CALIBRATED  OBLIQUE  I ON OG RAM  18. 
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Fig.  22.  EXPLP I MENTAL  CALIBRATED  OBLIQUE  IONOGRAM  22 
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EXPERIMENTAL  CALIBRATED  OBLIQUE  IONOGRAM  23. 


233 


Ill  POSSIBLE  MODES 


Several  possible  combinations  of  reflections  are  possible.  Two  of 

the  simplest  are  the  E-F  or  "n"  mode  and  the  FEF  or  "M"  mode,  as  shown 

in  Figs.  25  and  26.  (Henceforth,  within  the  multi-mode  terminology,  the 

sporad  j.c-E-layer  will  be  denoted  by  E,  rather  than  by  E  .  )  The  terminol- 

s 

ogy  adopted  for  the  description  of  the  modes  is  consistent  with  that 
adopted  at  the  meeting  in  Lindau,  May  1963  [Ref.  16,  p.  189],  and  uses  a 
dash  to  indicate  the  occurrence  of  a  ground  reflection.  Thus,  for 
example,  the  notation  E-F  (see  Fig.  25)  denotes  first  a  reflection  from 
the  Eg-layer,  then  a  reflection  from  the  ground,  followed  by  a  reflection 
from  the  F-layer  in  either  direction.  This  path  resembles  a  2-hop  F-mode 
path,  but  because  of  the  lower  height  of  the  Eg-layer,  the  time  delays 
for  the  E-F  mode  must  be  less  than  those  for  the  2-hop  F-mode;  furthermore 
since  the  angles  of  incidence  into  the  F-layer,  0  are  greater  than 
those  for  the  2F  case,  the  maximum  frequency  that  supports  the  E-F  mode 
(MUF)  should  be  higher  than  that  for  the  2F  mode.  These  expectations 
were  supported  by  the  experimental  records  obtained  . 

Similarly,  as  shown  in  Fig.  26,  the  FEF  mode  denotes  a  reflection 
from  the  F-layer,  followed  by  a  reflection  from  the  top  of  the  E  -layer, 
again  followed  by  a  reflection  from  the  F-layer,  which  directs  the  wave 
down  to  the  receiver. 

Other  possible  multi-modes  are  (refer  to  Figs.  27-32):  2E-F,  2F-E, 

FEFEF,  3E-F,  2E-2F,  FEF-E  (which  is  identical  to  2F  if  no  tilts  are 
present),  and  many  others  that  neither  are  very  likely  (e.g.,  EFEFE),  nor 
were  visible  within  the  time  delays  recorded  on  the  experimental  records 
(e.g.,  3F-E) .  Occurrence  of  the  FEF-F  mode  is  also  likely;  however,  it 
was  not  included  in  the  present  analysis. 
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Fig.  25.  MULTI-MODE  DIAGRAM:  E-F . 


Fig.  26.  MULT I -MODE  DIAGRAM:  FEF . 
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Fig.  27.  MULTI-MODE  DIAGRAM:  2E-F . 
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IV  THE  IONOSPHERIC  MODELS 


A .  The  Sporadic-E-Layer 

Evidence  of  a  sporadic-E-layer  is  shown  by  the  horizontal  traces  on 
the  records  from  1130  hr  on;  sane  of  the  records  in  Figs.  1-24  are  labeled 
accordingly.  The  evidence  is  even  stronger  on  the  records  of  1830  on. 

A  trace  with  this  behavior  is  caused  by  an  abrupt,  very  dense  layer  of 
electrons;  this  type  of  ionospheric  layer  permits  very  little  ray  pene¬ 
tration  into  the  layer  as  the  frequency  is  increased,  and  it  usually 
enables  the  achievement  of  a  rather  high  MUF . 

Rays  are  permitted  to  pass  to  the  F-layer  over  a  wide  range  of 
take-off  angles,  as  seen  for  1-,  2-,  and  3-hop  F  traces.  Clearly, 
either  the  E^layer  is  confined  to  the  center  of  the  path  and  is  limited 
in  extent,  or  it  extends  perhaps  over  the  entire  path,  but  consists  of 
irregularly  sized  and  spaced,  highly  dense  blobs  of  electrons.  The 
latter  would  permit  some  of  the  wavefront  energy  to  pass  through  the 
layer,  while  causing  some  of  it  to  be  reflected  by  the  highly  dense 
blobs.  The  author  is  inclined  to  accept  this  description  of  the  layer, 
especially  at  the  later  times,  since  the  presence  of  many  different 
multi-modes  suggests  a  fairly  large  E^  patch,  capable  of  supporting 
reflections  all  along  the  path.  Furthermore,  it  was  also  clear  that  2E 

s 

was  present,  suggesting  a  horizontal  layer  extent  of  at  least  700  km. 

This  size  is  consistent  with  other  findings  (e.g.,  see  Ref.  21). 

The  thickness  of  the  E  -layer  may  be  on  the  order  of  one  to  several 
kilometers.  J.  Carl  Seddon  [Ref.  22,  p.  78]  conducted  a  study  that 
utilized  measurements  of  electron  density  profile  made  by  comparing  the 
relative  doppler  shift  of  harmonically  related  signals  radiated  by  a 
vertically  ascending  rocket.  T.  Croft  [Ref.  23]  has  also  compiled  a 
set  of  records  containing  rocket  measurements  of  ionospheric  density, 
which  involve  E^-layers.  Both  of  these  measurements  confirm  that  the 
average  "thickness"  of  the  layer  is  about  2  km,  and  heights  in  the 
vicinity  of  100  km  occur.  Typical  profiles  taken  from  Ref.  22  are 
reproduced  in  Fig.  33.  The  layer  thickness  is  measured  between  points 
where  the  density  is  several  orders  of  magnitude  below  the  maximum  den¬ 
sity.  Since  a  wave  will  penetrate  the  E^layer  to  a  small  degree — 
either  from  above  or  from  below — the  effective  thickness  of  the  layer 
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(a)  Comparison  of  sporadic  E  at  noon  over  New  Mexico, 

U.S.A  with  that  over  Woomera,  Australia  at  midnight. 


ELECTRON  DENSITY  (I03el/ee) 

(b)  Detailed  profile  of  sporadic  E  over  Ft.  Churchill. 

Fig.  33.  ELECTRON  DENSITY  PROFILES  MEASURED  WITH  ROCKETS.  (From  Smith 
and  Matsushita,  Ionospheric  Sporadic  E,  The  Macmillan  Co.,  New  York, 
1962 .  Page  78. ) 
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is  probably  on  the  order  of  1  km.  However,  since  the  path  length  was 
appreciably  greater  than  this  thickness,  it  was  assumed  for  greater 
simplicity  that  the  layer  was  of  zero  thickness. 

We  then  arrive  at  the  following  preliminary  model  assumptions: 

(1)  large  extent,  possibly  over  the  entire  path; 

(2)  very  thin,  highly  dense  blobby  layer; 

(3)  nonblanketing  (holes  between  blobs). 

The  height  of  the  Eg-layer  was  determined  through  knowledge  of  the 
time  delay  of  propagation  over  the  path,  and  by  assuming  curved-earth 
geometry.  The  time  delay,  denoted  hereafter  by  T  ,  and  height  are  relat¬ 
ed  by 


D 


(for  1-hop)  , 


(for  2-hop)  , 


where 

R  =  mean  radius  of  earth,  6370  km, 

Dq  =  path  length,  1881.27  km, 

h  =  height  of  E  -layer. 

s 

These  functions  were  plotted  for  this  particular  path,  as  shown  in  Fig. 

34.  The  T  's  corresponding  to  the  IE  and  2E  modes  were  extracted  from 
u  s  s 

the  records  and  are  tabulated  in  Table  1 .  The  corresponding  E  heights 
vs  midpath  time  of  day  are  plotted  in  Fig.  35;  there  are  two  curves,  one 
for  1-hop,  the  other  for  2-hop.  For  the  later  times,  the  2-hop  data 
indicated  lower  heights  than  those  found  in  the  1-hop  data.  Ordinarily, 
one  would  expect  the  2-hop  mode  to  penetrate  slightly  deeper  into  the 
layer,  and  thus  be  reflected  from  a  higher  "virtual  height"  than  the 
1-hop  mode — as  seemed  to  be  the  case  earlier  in  the  day.  In  most  cases, 
the  2-hop  trace  appeared  at  lower  frequencies,  and  the  1-hop  at  higher 
frequencies;  at  these  higher  frequencies,  there  may  have  been  deeper 
1-hop  mode-layer  penetration.  However,  a  more  reasonable  explanation 
for  the  inconsistency  might  be  the  existence  of  two  or  more  different 
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levels  of  E  ionization.  Such  an  occurrence  has  been  reported  by  observ- 
s 

ers  in  the  past,  and  again  just  recently  at  Stanford  by  L.  E.  Sweeney, 

who  is  currently  taking  oblique-sounding  records  between  Los  Banos, 

California,  and  Bearden,  Arkansas,  using  a  very  large  aperture  receiving 

array.  Seddon  [Ref.  22,  p.  83]  cites  examples  in  which  "preferred" 

levels  of  ionization  separated  by  6  km  have  been  found.  The  existence 

of  these  levels  correlated  with  the  occurrence  of  wind  shear  levels 

separated  by  the  same  6  km.  It  is  therefore  possible  that  a  lower  level 

was  formed  over  midpath  later  in  the  day,  accompanied  by  the  disappearance 

of  levels  over  other  parts  of  the  path. 

From  Fig.  35,  it  is  impossible  to  tell  whether  an  E^-layer  of  higher 

(or  lower)  density  occurred  at  the  transmitter,  middle,  or  receiver  end 

of  the  path.  Therefore,  given  a  specific  ray  path  for  a  particular  mode 

from  which  an  E  reflection  occurs  at  a  specific  range,  the  exact  height 
s 

of  the  E  -layer  cannot  be  deduced.  An  average  E  height  was  therefore 
s  s 

determined  by  averaging  the  1-  and  2-hop  E  data  of  Fig.  35  at  each  time 

s 


MIDPATH  LOCAL  TIME 

Fig.  35.  SPORADIC-E  HEIGHT  VS  MIDPATH  LOCAL  TIME  (MLT ) . 
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of  day;  this  is  shown  in  Fig.  36,  Furthermore,  to  make  the  models  more 
consistent,  all  these  values  were  averaged  over  the  day's  data,  and  an 
average  heigut  of  117  km  was  obtained;  this  height,  by  coincidence,  is 
equal  to  one  of  Seddon's  "preferred"  heights.  The  height  of  117  km 
proved  quite  successful  for  model  matching,  except  at  1845  MLT,  in  which 
case  a  height  of  105  km  gave  a  much  better  fit.  This  exception  could 
have  been  caused  by  the  occurrence  of  a  small,  lower  E  patch,  located 
over  the  appropriate  reflection  points  along  the  path. 


Fig.  36.  AVERAGE  HEIGHTS  OF  FIGURE  35  VS  MIDPATH  LOCAL  TIME. 
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B. 


The  F-Layer  Models 


1 .  The  Initial  Computer  Models 

One  of  the  simplest  possible,  yet  realistic,  models  for  the 
F-layer  is  the  (3-Chapman  layer,  given  by 

h  -  h 

N  =  N  exp(l  -  z  -  e  )  ,  where  z  =  — — -  , 

o  H 

s 

and 

N  =  electron  density, 

Nq  =  maximum  electron  density, 
h  =  height, 

h  =  height  of  maximum  ionization, 
m 

H  =  scale  height  in  (3-Chapman  equation  for  density, 
s 

This  model  was  quite  successful,  and  with  suitable  modifications  in  N  , 

o 

h^,  and  Hg,  was  used  throughout  the  analysis.  Table  2  lists  the  parame¬ 
ters  used  in  the  generation  of  these  initial  computer  models. 


Table  2 


F-LAYER  IONOSPHERES  TRIED  IN  THE  INITIAL  MATCHING  PROCEDURE 


Model 

Maximum 

Density  N 

o 

(electrons/cm^ ) 

Max imum 
Height  h 

m 

(km) 

Sc  a  lo 

Height  H 

s 

(km) 

A 

io6 

270 

30 

B 

io6 

270 

40 

C 

106 

300 

40 

D 

106 

300 

65 

E 

106 

300 

90 

F 

io6 

350 

40 

G 

106 

350 

65 
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2. 


Modification  of  the  Models  To  Match  the  MUFs 


The  synthesized  ionograms  are  presented  as  plots  of  T^  vs 
frequency .  The  parameter  that  governs  the  behavior  of  a  wave  propagating 
through  the  ionosphere  is  the  index  of  refraction,  given  by 


2 

M  = 


1 


80.6 


N  (h) 


N  g (h) 

=  1  -  80.6  -£ -  , 

f2 

where 

N (h)  4  N  g (h)  , 
o 

g(h)  =  profile  shape  factor, 
f  =  frequency  in  MHz. 

t 

Thus,  if  the  density  factor,  N  ,  is  changed  to  N  =  kN  ,  where 
2  ^  o  o 

k  is  some  constant,  then  ^  will  be  left  unchanged,  providing  the 

frequency  is  scaled  by  the  factor  Vk\  Correspondingly,  if  we  are  given 
a  synthesized  ionogram  on  which  the  frequency  scale  is  changed  by  the 
factor  \lk,  the  equivalent  model  density  is  changed  again  by  the  factor 
k.  Therefore,  to  match  the  records,  the  model's  frequency  was  scaled  to 
give  an  MUF  identical  to  the  experimental  MUF;  both  cases  are  for  the  IF 
mode.  It  is  perhaps  easiest  to  visualize  this  by  plotting  the  ionogram 
on  a  log' frequency  scale;  the  scale  may  then  be  shifted  along  the  hori¬ 
zontal  axis  to  yield  an  ionogram  for  different  values  of  N  .  Since  the 

o 

experimental  data  were  recorded  on  linear  scales,  the  shifting  procedure 
was  not  feasible;  no  attempt  was  made  to  plot  the  experimental  data  on 
semilog  coordinates. 
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V  COMPUTER  SYNTHESIS 


The  models  for  the  sporadic-E-  and  F-layers  mentioned  above  were 
employed  in  the  generation  of  computer  ionograms .  The  first  task  was  to 
find  a  model  that  yielded  an  ionogram  whose  IF  mode  matched  the  corre¬ 
sponding  IF  mode  of  an  experimental  record. 

Curved-earth  geometry  was  employed,  and  horizontal  tilts  in  electron 
density  were  ignored.  The  sporadic-E-layer  was  assumed  to  extend,  over 
the  entire  path  at  a  constant  height.  ad  was  made  to  be  partially 
reflecting  and  partially  transmitting  ;  ,vith  time  retardation  ignored). 

The  ray-tracing  procedure  (denoted  Ray-319)  was  developed  by 
T.  Croft  and  P.  Fialer  at  the  Stanford  Radioscience  Laboratory.  Given 
the  N(h)  profile,  the  program  specifies  a  frequency  and  take-off  angle, 
and  then  follows  the  ray  through  the  ionosphere.  The  refractive  index 
is  calculated  at  1-km  intervals;  then  Snell's  law  is  solved  at  each 
successive  level.  The  ray  is  followed  until  it  strikes  the  ground,  at 
which  point  its  horizontal  distance  from  the  transmitter  is  recorded; 
the  group  and  phase  time  delays  and  the  true  height  of  reflection  are 
also  recorded.  A  new  take-off  angle  is  then  specified,  and  the  process 
is  repeated.  The  output  was  recorded  on  the  line  printer  and  was  also 
put  on  punched  cards  for  later  use  in  an  interpolation  program. 

A .  Calculations 

Given  a  particular  mode,  the  height  of  the  Eg-layer,  a  model  F-layer, 
and  a  particular  frequency,  there  are  only  two  possible  take-off  angles 
that  will  support  the  propagation  of  a  ray  over  the  fixed  distance 
between  the  transmitter  and  receiver,  namely,  the  lower  and  upper 
(Pederson)  rays  in  the  F-layer  (magnetoionic  splitting  is  ignored). 

Let  us  concentrate  on  the  lower  ray.  To  illustrate  the  method  of 
calculation,  we  will  consider  the  geometry  presented  by  the  E-F  or  "N" 
mode,  shown  in  Fig.  37  .  As  seen,  the  F-layer  has  been  replaced  by  its 
equivalent  virtual  height  of  reflection,  h’,  which  corresponds  to  the 
take-off  angle  solution,  p,  for  the  lower  ray;  h'  will  be  a  function  of 
the  frequency  and  ionospheric  properties  of  the  F-layer.  We  then  see 
that,  for  a  given  E-layer  height,  h*  there  will  be  a  specific  distance, 
Dg,  which  exists  for  the  solution  p.  Moreover,  there  will  be  a  specific 
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Fig.  37.  GEOMETRY  OF  E-F  MODE  FOR  CALCULATION  OF  TIME  DELAYS. 


distance,  D  ,  over  which  the  F  ray  propagates.  For  the  E-F  mode, 
F 


D  +  D  =  D  . 
E  F  o 


Dq  was  determined  through  knowledge  of  the  latitude  and 


longitude  of  the  transmitter  and  receiver,  and  found  to  be  1881.27  km. 
D£  is  a  function  of  p,  and  is  given  by 


de  =  2R 


h'1  (^)  ■ p] ■ 


where 

R  =  6370  km  as  before, 

h'  =  specified  117  km  in  all  cases  but  one,  where  105  km  was  used 
E 

Thus,  the  remaining  distance  D  is  given  by  D =  D  -  D  .  Corresponding 

r  r  O  F 

ly,  the  ray-tracing  program  presents  an  output  of  skip  distance  vs  (3. 

We  shall  call  this  distance  D^;  the  solution  p  exists  where  D^  =  D^ . 
Thus,  D  was  tabulated  for  values  of  p  ascending  in  short  steps;  then, 
both  this  table  and  the  computer  output  were  visually  scanned  for  the 
approximate  solution  p,  which  was  flanked  by  two  values  of  p,  above 
and  below,  in  each  table.  The  cards  representing  these  values  were 
inserted  into  a  program  that  interpolated  for  the  correct  solution  p; 
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this  process  is  shown  graphically  in  Fig.  38,  where  the  two  distance 

curves  are  plotted  vs  take-off  angle  and  solutions  are  found  for 

different  frequencies.  The  solution  p  was  then  used  to  find  the  solution 

time  delay  for  the  F  hop;  finally  the  total  time  delay  over  the  patn  at 

each  particular  frequency  was  determined.  For  the  E-F  case  T  =  T  j-  T 

’  D  F  E  ‘ 

At  those  frequencies  for  which  an  upper  ray  was  supported,  the  process 
was  repeated  at  a  higher  set  of  p's  (Fig,  38).  The  frequencies  were 
chosen  1  MHz  apart;  the  set  of  solutions  obtained  for  all  the  frequencies 
was  plotted,  and  the  points  were  connected  to  form  the  ionogram  trace. 

The  solution  was  always  a  smooth  curve,  since  the  F-layer  profile  was 
a  monotonic  function. 


Fig.  <8.  F-LAYER  GROUND  DISTANCE  VS  TAKE-OFF  ANGLES, 
I L'  STRATI NG  COMPUTER  INTERPOLATION, 
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Before  we  generalize  the  procedure,  consider  for  a  second  example 
the  FEF,  or  "m"  mode,  shown  in  Fig.  39.  For  this  mode, 


°F  = 


D  +  D 
o  E 


The  function  D  is,  of  course,  the  same  as  before.  However,  a  new  set 

of  curves  D^(p)  is  obtained.  The  graphical  (computer-interpolation) 

procedure  is  carried  out  as  before.  Finally,  the  time  delay  is  obtained 

as  T  =  2T  -  T  .  Note  that  the  E  hop  is  fictitious;  yet  it  is  included 
D  F  E 

to  complete  the  geometry. 


Dp  3  (D0  +  De)/2 

tD  *  2TF  -  te 

Fig.  39.  GEOMETRY  OF  FEF  MODE  FOR  CALCULATION  OF  TIME  DELAYS. 

The  same  procedure  may  be  repeated  for  all  other  possible  modes. 

For  each  mode,  we  will  have  a  new  set  o*  functions  for  D  (p),  calculated 

from  knowledge  of  the  geometry.  These  new  functions  will  intersect  the 

computer  ray  tracings  D^(p),  as  before,  and  a  new  ionogram  trace  is 

finally  obtained  by  connecting  the  time-delay  points  obtained  at  each 

new  frequency.  Table  3  lists  the  various  possible  modes,  together  with 

their  functions  for  D  and  T  .  Note  that  some  modes  could  be  expressed 

F  D 
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Table  3 


differently;  for  instance,  2E-2F  could  be  2F-2E,  E-F-E-F,  F-E-F-E, 

E-2F-E,  or  F-2E-F .  If  tilts  in  electron  density  were  present,  these 
modes  would  be  separated;  otherwise,  they  would  all  appear  superimposed. 

B .  Matching  the  Models  to  IF  and  2F  Traces 

The  model-matching  procedure  is  essentially  a  "cut-and-try "  process; 

however,  certain  guidelines  may  be  followed  in  order  to  make  a  good  guess 

at  the  first  "cut."  The  following  is  an  account  of  how  the  models  used 

were  derived  through  trying  different  values  of  h  ,  H  ,  and  the  modified 

m  s 

N  . 
o 

Model  E  (Table  2)  was  tried  first,  since  it  had  been  used  on  numer¬ 
ous  occasions  in  other  studies  by  T.  Croft  at  Stanford  (the  Stanford 
number  is  I  ID  165).  The  TQ  at  the  MUF  of  model  E  closely  matched  the 
MUF  of  Ionogram  14  (Fig.  14);  however,  the  MUFs  were  different.  There¬ 
fore,  the  frequency  scale  of  the  model  was  modified  according  to  the 
technique  mentioned  above  : 


Model  MUF  =  20.8  MHz. 
Experimental  MUF  =13.8  MHz. 


Therefore,  Vk  6  =  13.8/20.8  =  0.663.  The  frequency  scale  of  model  E 

was  multiplied  by  this  factor.  The  result  is  shown  in  Fig.  40.  As  seen, 

at  lower  frequencies  the  model's  time  delays  are  too  low.  The  scale 

height  was  reduced  to  H  =40  km,  and  the  maximum  height  was  increased 

s 

to  h  =  350;  model  F  resulted.  The  time  delay  at  the  MUF  was  7.0  msec, 
m 

which  was  much  too  high,  even  for  records  at  later  times  of  the  day. 

The  scale  height  was  increased  to  65,  and  model  G  resulted;  unfortunately, 
the  MUF  T^  was  still  too  high. 

It  was  noted,  however,  that  when  the  scale  height  is  decreased,  the 

ionization  becomes  more  concentrated  at  heights  near  h  ,  which  raises  the 

m 

T  's  at  the  lower  frequencies.  Therefore,  h  was  returned  to  300  km,  and 
D  m 

H  was  left  at  65  km;  model  D  resulted.  When  this  model's  MUF  was  matched 
s 

to  Ionogram  14,  a  successful  match  was  obtained.  was  found  to  be 

0.640,  and  the  new  model  is  given  the  number  Dl,  where  it  differs  from 

2  5 

model  D  by  its  maximum  density,  IT  =  Nq  =  4.09(10)  . 

It  was  also  noted  that  Ionogram  15  gave  time  delays  somewhat  higher 

at  lower  frequencies;  however,  the  Tq  at  its  IF  MUF  was  approximately  the 

same  as  that  of  Ionogram  14.  Hence,  hm  was  left  at  300  km,  and  Hg  was 

lowered  to  40  km;  model  C  resulted,  whose  MUF  was  scaled  to  Ionogram  15 

2  5 

with  K.  =  0.622.  This  yielded  model  Cl,  with  its  N’  =  K  N  =  3.87(10)  . 
f  o  f  o 

Third,  note  that  the  MUF  T^'s  are  somewhat  lower  at  earlier  times  of 

the  day,  as  in  Ionograms  9  through  13.  Therefore,  since  it  was  apparent 

that  the  drop  in  MUF  T  could  be  realized  by  lowering  h  ,  the  height  of 

D  m 

maximum  ionization  was  dropped  to  270  km.  Two  values  of  H  ,  30  and  40  km, 

s 

were  tried,  and  models  A  and  B  resulted.  Model  A  yielded  too  high  time 
delays  at  lower  frequencies.  Model  B,  however,  through  proper  adjustment 
of  the  frequency  scale,  gave  good  matches  to  the  IF  traces  of  Ionograms  9, 
12,  and  13.  The  results  were 


Ionogram 

9,  Kf 

=  0.816, 

N  ’  =  6.65(10)  — 
o 

-model 

B1 

Ionogram 

12,  Kf 

=  0.705, 

N'  =  4 . 96 (10 ) 5  — 
o 

-model 

B2 

Ionogram 

13,  Kf 

=  0.654, 

N'  =  4,27(10 )5 — 
o 

-model 

B3 
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The  five  models,  Bl,  B2,  B3,  Dl,  and  Cl,  are  plotted  in  Figs.  41  and  42. 
Note  the  change  in  the  F-layer  ionosphere  as  the  day  progresses . 

C.  Including  the  E_-Layer:  Generation  of  Multi-Modes 

“ - S - —  1  ■  —  - ■  -  . . . . . — 

For  lonograms  12  through  15,  an  Eg-layer  height  of  117  km  was 
employed  in  the  programs  used  to  generate  the  multi-modes.  For  Ionogram 
9,  an  E^  height  of  105  km  gave  a  much  closer  match.  The  effect  of 
varying  this  height  may  be  seen  in  Chapter  VII,  Error  Analysis.  The 
results  are  shown  in  Figs.  43-47.  They  are  discussed  individually  below. 


Fig.  41.  PROFILES  OF  IONOSPHERIC  MODELS  USED:  MODELS  Bl,  B2,  AND  B3 . 
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Fig.  42.  PROFILES  OF  IONOSPHERIC  MODELS  USED:  MODELS  D1  AND  Cl. 
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RESULTS  OF  IONOGRAM  MATCHES:  I ON OG RAM  9,  MODEL  B1 


RESULTS  OF  IONOGRAM  MATCHES:  IONOGRAM  12,  MODEL  B2 


RESULTS  OF  IONOGRAM  MATCHES:  IONOGKAM  14,  MODEL  D1 
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RESULTS  OF  IONOGRAM  MATCHES:  IONOGRAM  15,  MODEL  Cl 


D  . 


Observations 


1 •  lonogram  9;  1845  MLT  (Fig.  43) 

This  was  the  best  match  obtained.  As  seen,  modes  IF,  E-F,  FEF 

2F-E,  and  2F-2E  are  clearly  explained.  However,  it  is  not  clear  what 

occurred  near  the  2F  mode.  For  frequencies  below  about  11  MHz,  it 

appears  that  a  2F  or  FEF-E  mode  is  present.  At  11  MHz,  a  peculiar  split 

ting  of  the  mode  suggests  that  both  2F  and  FEF-E  are  present;  2F  is  the 

most  probable,  since  we  are  not  sure  about  the  extent  of  the  E  -layer. 

Therefore,  che  trace  at  the  lower  frequencies  is  probably  2F,  while  the 

trace  above  11  MHz  is  probably  FEF-E.  Furthermore,  note  that,  if  the 

FEF  portion  of  the  latter  mode  occurred  at  the  east  end  of  the  path,  its 

higher  time  delays  would  be  explained  because  the  electron  density  for 

the  two  F  hops  would  be  slightly  lower  than  for  the  2F  mode.  (The 

subject  of  tilts  is  discussed  later.) 

As  noted,  an  E  height  of  105  km  was  employed  in  this  model, 
s 

Reference  to  Fig.  35  does  not  suggest  this  height.  However,  when  117  km 

was  used,  there  was  a  fairly  large  discrepancy  between  the  model  and  the 

record  trace.  As  mentioned  in  Chapter  IV-A,  a  small  E  patch  at  a  lower 

s 

height  may  have  occurred  somewhere  along  the  path;  or  an  error  in  time- 
delay  calibration  on  the  experimental  records  caused  an  error  in  the 
estimation  of  E^  virtual  heights. 

2 .  lonogram  12:  1945  MLT  (Fig.  44) 

This  was  also  a  fairly  good  match.  The  modes  IF,  E-F,  FEF, 
and,  seemingly,  2F-E  are  explained.  The  occurrence  and  splitting  of  the 
modes  2F  and  FEF-E  are  also  apparent. 

The  group  of  modes  around  2F-2E  could  very  well  be  modes  such 
as  F-E-F-E  and  2E-2F  (as  mentioned  in  Chapter  III)  which  have  the  sane 
model,  but  which  see  slightly  different  ionospheres  in  ray  travel. 
Because  it  js  even  possible  that  3E-F  was  present  in  this  group,  this 
mode  is  dotted-in  for  reference;  however,  it  is  not  probable,  since  3E 

s 

was  not  evident  on  the  records. 
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3. 


lonogram  13;  2000  MLT  (Fig.  45) 


This  was  not  a  very  exceptional  match;  nevertheless,  we  may 

still  identify  the  IF,  2F,  E-F,  and  FEF  modes.  It  seems  also  clear  that 

the  2F-E  mode  was  present.  We  note  that  IE  and  2E  are  becoming  more 

s  s 

pronounced  in  this  and  the  following  records. 

4.  lonogram  14:  2030  MLT  (Fig.  46) 

Modes  IF,  E-F,  and  FEF  are  clearly  identifiable.  Modes  2F 
and/or  FEF-E  also  appear  with  2F  again  the  most  probable. 

5 .  lonogram  15:  2100  MLT  (Fig.  47) 

Although  this  was  a  poor  match,  we  may  still  label  the  IF,  E-F, 
and  FEF  modes.  Because  of  the  model/experiment  discrepancy  for  IF,  we 
may  identify  2F  as  the  major  trace  just  above  the  FEF  mode,  since  the 
model  gives  time  delays  too  high  at  the  lower  frequencies. 

Note  some  general  characteristics  that  occur  on  each  model  ionogram 
with  respect  to 

(1)  the  general  location  of  the  modes  with  respect  to  IF,  2F,  and 
3F  traces  (see  Chapter  VI); 

(2)  the  extent  of  the  upper  ray  into  higher  time  delays,  which  is 
determined  by  the  distribution  of  ionization  near  the  maximum 
density; 

(3)  the  fact  that  if  the  upper  ray  of  the  IF  trace  were  extended  to 

higher  time  delays,  to  the  top  of  the  record,  and  if  no  tilts 

were  present,  all  mode  traces  would  appear  to  the  left,  and 

above,  the  IF  trace — except  the  IE  and  2E  modes,  which  do 

s  s 

not  depend  on  the  F-layer. 
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VI  COMPARISON  WITH  SIMPLER  MODEL 


To  measure  the  effectiveness  of  the  aforementioned  model-comparison 
technique,  a  simpler  model  was  used  to  predict  the  FEF  and  E-F  mode 
structure.  Both  the  E  -  and  F-layers  were  assumed  to  be  copper  sheets-- 
the  mirror  analogy;  additionally,  the  Eg-layer  was  allowed  to  transmit 
part  of  the  energy  incident  upon  it,  so  that  it  was  nonblanketing. 

The  E^-layer  was  assumed  to  be  117  km  high;  the  height  of  the 
F-layer  was  varied.  Curved-earth  geometry  was  assumed,  and  the  appropri¬ 
ate  formulas  were  programmed  into  a  digital  computer  to  calculate  the 
time  delay  and  F-layer  virtual  height  which  would  support  propagation  of 
a  particular  mode  at  a  given  take-off  angle.  A  similar  analysis  was 
conducted  by  Bowhill  [Ref.  24]  in  1961.  The  formulas  used  are  as  follows 

The  take-off  angle  (3  is  specified;  then 

De  =  2R  ^ cos*"1  (f~^j  ~  p]  ^  before, 

where 

R  =  6370  km, 
hi  =  117  km. 

hi 

Also 


=  <P, 


X  Do  +  P2  X  V/P3  ' 


where  for  each  mode,  P  ,  P  ,  and  P  are  defined  by  the  corresponding 

1  ^  u 

equation  for  D  appearing  in  Table  3. 

F 

Then 


msec  , 


msec , 


Tn  =  Q-I  XTr  K!„X  T 


msec, 
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TIME  DELAY  (msec) 


where  Q  and  Q  are  specified  for  each  mode  according  to  the  T  equations 

1  Z  L) 

in  Table  3.  The  time  delays  vs  F-layer  virtual  height  were  plotted  as 
shown  in  Fig.  48. 

To  predict  the  time  delays  at  which  a  particular  mode  would  occur, 
it  was  necessary  first  to  choose  a  frequency,  and  then  to  estimate  the 
average  virtual  height  for  the  F-layer,  from  the  IF  and  2F  modes,  accord¬ 
ing  to  Fig.  48.  We  note  that  the  2F  virtual  height  will  be  greater  than 


Fig.  48.  TIME  DELAY  VS  F-LAYER  VIRTUAL  HEIGHT  FOR  SEVERAL  MODES— 
MIRROR-ANALOGY  MODEL. 
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the  IF  virtual  height,  since  in  the  latter  case  the  wave  does  not 
approach  the  layer  from  as  steep  an  angle.  Table  4  lists  these  results. 


Table  4 

IF  AND  2F  VIRTUAL  HEIGHTS  (h^)  FOR  EACH  IONOGRAM  AT  SPECIFIC  FREQUENCIES 


Ionogram 

Frequency 

(MHz) 

V1F 

(msec ) 

Td-2F 

(msec ) 

hp-lF 

(km) 

h;-2F 

(km) 

hp  (avg) 

(km) 

9 

9 

6.57 

7.13 

237 

238 

238 

12 

9 

6 .58 

7.23 

240 

253 

247 

13 

8 

6  .56 

7.25 

232 

257 

245 

14 

7 

6.54 

7.30 

221 

264 

243 

15 

7 

6.60 

7.44 

250 

284 

267 

The  F-layer  virtual  heights  were  then  used  to  predict  time  delays 
for  the  FEF  and  E-F  modes.  Also,  at  the  frequencies  specified  for  each 
ionogram,  the  time  delays  for  the  two  modes  (FEF  and  E-F)  appearing 
between  the  IF  and  2F  modes  were  recorded.  The  results  appear  in  Table  5  . 


Table  5 

COMPARISON  OF  MIRROR -ANALOGY  MODEL  AND  EXPERIMENTAL  IONOGRAM  TIME  DELAYS 
AT  SPECIFIC  FREQUENCIES  OF  TABLE  4,  USING  h^ (average) 


Ionogram 

Frequency 

(MHz) 

td  e“f 

(msec ) 

Tq  FEF 

(msec ) 

Calculated 

Observed 

Calculated 

Observed 

9 

9 

6.79 

6.77 

6.84 

6.85 

12 

9 

6.82 

6.78 

6.89 

6.87 

13 

8 

6.84 

6.78 

6.88 

6.87 

14 

7 

6.80 

6.80 

6.87 

6.91 

15 

7 

6.87 

6.80 

7.00 

7.00 
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It  is  clear  from  the  aforementioned  model  predictions  that,  at  the 

frequencies  considered,  the  lower  of  the  time  delays  observed  correspond 

to  the  E-F  mode.  This  fact  is  also  predicted  by  the  simpler  model,  as 

seen  in  Table  5 .  It  is  probably  most  reasonable  to  consider  the  average 

calculated  time  delays,  since  the  virtual  height  of  the  F-layer  that 

supports  propagation  of  the  E-F  and  FEF  modes  must  lie  somewhere  between 

the  IF  and  2F  virtual  heights.  This  may  be  easily  seen  by  noting  that 

the  equivalent  F-hop  distance,  D  ,  is  always  greater  than  |D  ,  but  less 

F  o 

han  D^,  for  both  modes,  for  all  virtual  heights  above  117  km.  (A 

longer  D  results  in  a  lower  virtual  height.) 

F 

All  but  Ionogram  14  have  given  FEF  predictions  within  0.02  msec, 
and  Ionogram  14  gives  it  within  0.04  msec.  The  predictions  for  the  E-F 
mode  are  almost  as  accurate:  within  less  than  0.04  msec  for  Ionograms 
9,  12,  and  14,  less  than  0.06  msec  for  Ionogram  13,  and  less  than  0.07 
msec  for  Ionogram  15. 

This  comparison  gives  further  justification  to  the  computer-modeling 
procedure,  and  could  even  be  used  as  a  rough  estimation  of  mode  structure, 
when  only  a  good  guess  at  the  type  of  mode  is  needed. 

The  presentation  of  Fig.  48  looks  quite  similar  to  an  oblique 
ionogram,  with  the  frequency  scale  replaced  by  F-layer  virtual  height. 

The  similarity  lies  in  the  fact  that  for  higher  frequencies,  the  virtual 
height  of  reflection  does  go  up,  because  refraction  by  the  ionosphere  is 
less  pronounced  at  these  higher  frequencies.  The  exact  relationship 
between  virtual  height  and  frequency  will  depend  on  the  electron  profile 
over  the  path  (or  the  model  assumed). 

Figure  48,  probably  one  of  the  most  significant  in  the  report,  shows 
in  a  general  way  how  modes  are  likely  to  be  placed  in  relationship  to 
each  other,  especially  when  a  single  F-layer  is  present.  The  similarity 
may  be  seen  through  comparison  with  the  model-matches  done  previously. 
Note  in  particular  the  FEF  and  E-F  modes  of  Fig.  46.  The  modes  cross  at 
about  5.3  MHz,  with  the  E-F  mode  remaining  the  lower  after  the  cross¬ 
over.  Exactly  the  same  effect  occurs  in  Fig.  48,  where  the  cross-over 

occurs  at  h'  =  220  km.  The  time  delay  at  the  cross-over  is  the  same  for 
F 

both  models — 6.76  msec.  Reference  to  the  models  in  Figs.  44  and  45 
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shows  the  same  cross-over.  In  Figs.  43  and  47,  the  cross-over  occurs  at 
a  much  lower  frequency,  if  at  all;  the  E-F  mode  is  still  the  lower  of  the 
two.  Precisely  the  same  effect  is  noted  for  the  2E-F,  2F,  and  FEFEF 
combination  in  Figs.  43  and  44,  when  compared  again  with  Fig.  48. 

The  simple  model  predictions  of  Fig.  48 — or  curve  sets  like  it  for 
different  ranges  and  E^  heights--should  provide  a  valuable  tool  for 
deducing  the  types  and  positions  of  modes.  On  a  long  path,  and  especial¬ 
ly  if  severe  tilts  occur,  the  simple  model  could  prove  to  be  entirely 
inadequate;  however,  a  relatively  good  "guess"  is  still  possible. 
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VII  ERROR  ANALYSIS 


Various  possible  errors  can  occur  within  the  mode-matching  or 
predicting  analysis.  Basically,  we  are  comparing  the  traces  on  the 
experimental  data  with  calculated  or  synthesized  traces  fro-a  a  model. 

The  following  problems  could  exist: 

(1)  Noisy  data:  ground  forward  scatter,  spread-F  (if  present), 
Faraday-rotation-with-f requency  blobs  (Ref.  25),  blanketing 
that  causes  cut-off  at  lower  frequencies,  excessive  interfer¬ 
ence  from  broadcast  stations,  and  other  interference. 

(2)  Time-delay  error  and  smearing  of  traces  due  to  propagation  of 
waves  through  a  blobby  E^-layer — poor  resolution. 

(3)  Improper  calibration  of  data  time-delay  scales. 

(4)  Misalignment  of  time-delay  scales  and  broadening  of  ticks  due 
to  photograph  scale  expansion — poor  resolution. 

(5)  Inadequate  model  for  F  and  E  layers.  Most  of  the  error  in 

s 

the  latter  is  attributable  to  errors  in  measurement  of  virtual 
height . 

The  errors  due  to  (1)  and  (2)  are  very  difficult  to  discuss  quanti¬ 
tatively.  Those  due  to  (3)  are  assumed  not  to  exist,  although  such  errors 
could  explain  why  a  lower-than-expected  Eg  height  gave  a  better  match 
in  Ion ogram  9,  and  would  probably  also  yield  closer  matches  in  the  other 
ionograms.  The  inaccuracies  caused  by  (4)  were  minimized  by  care  in 
data  reproduction  and  model-experiment  scale-matching.  The  errors  due 
to  (5)  are  probably  the  greatest,  and  may  be  discussed  quantitatively. 

The  F-layer  model  could  be  modified  to  any  desired  extent.  A 
parabolic  layer  might  be  tried;  it  has  given  good  matches.  Extra  elec¬ 
trons  could  be  substituted  into  any  model  on  a  trial-and-error  basis  to 
match  to  the  IF  and  2F  traces  more  closely .  This  is  a  very  slow  process, 
although  the  technique  could  be  improved  by  having  the  computer  do  all 
the  interpolation  procedures,  etc.,  so  that,  given  inputs  of  distance, 
Es~layer  height,  earth  radius,  F-layer  model,  and  type  of  mode  desired, 
an  ionogram  output  is  produced.  One  could  also  tell  the  computer  what 
the  desired  MUF  is,  and  it  would  modify  the  frequency  scale  accordingly. 

A  most  extravagant  idea  would  be  to  devise  a  self -adapting  program  to 
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match  to  any  desired  accuracy. 

To  estimate  the  error  in  time  delay  due  to  an  error  in  E^-layer 
height,  ionosphere  model  A  (Table  2)  was  used  in  the  calculations  of  IF, 
2F,  E-F,  FEF,  2E-F,  and  FEFEF  modes.  Heights  of  105,  117,  and  124  km 
were  used  to  produce  3  ionogram  sets,  as  shown  in  Figs.  49-51.  Compari¬ 
son  of  the  figures  is  self-explanatory.  An  example  of  the  uncertainty 
in  time  delay  may  be  calculated  at  11  MHz  for  the  FEF  mode,  where  the 
time  delay  for  E^-layer  heights  of  105,  117,  and  124  km  are  6.90,  6.87, 
and  6.84  msec,  respectively.  If  a  height  of  117  km  is  assumed  (as  it 
was),  then  an  error  of  0.03  msec  will  occur  if  the  height  is  actually 
105  or  124  km.  This  T^  error  is  of  the  same  order  obtained  with  the 
simplified  model,  and  is  also  of  similar  order  to  that  observed  in  the 

Chapman  model-matching.  Note  that  over  the  day  of  measurement,  the  E  - 

s 

layer  height  varied  between  107  and  124  km. 

There  are  also  some  approximation  errors  in  the  calculations  due  to 
the  linear  interpolation  done  in  the  computer;  however,  these  are  of  no 
great  significance  compared  with  the  resolution  obtained  in  the  data  of 
Figs.  1-24,  which  are  limited  by  factors  of  equipment  and  data  reproduc¬ 
tion  . 
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Fig.  49.  ERROR  ESTIMATE:  MODEL  IONOGRAMS  FOR  SPORADIC-E  HEIGHTS  OF  105 
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Fig.  50.  ERROR  ESTIMATE:  MODEL  IONOGRAMS  FOR  SPORADIC-E  HEIGHTS  OF  117 
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Fig.  51.  ERROR  ESTIMATE:  MODEL  IONOGRAMS  FOR  SPORADIC-E  HEIGHTS  OF  124  KM. 


VIII  TILTS  IN  ELECTRON  DENSITY:  MODELING  SUGGESTIONS 


It  is  highly  probable  that  the  electron  density  profile  was  not  at 
all  times  constant  over  the  E-W  path  between  Stanford  and  Lubbock,  since 
the  sun  s  zenith  angle  was  different  over  the  path, 

Because  each  mode  propagates  chrough  different  sections  of  the 
ionosphere,  a  simple  laterally  constant  model  would  not  be  adequate  to 
explain  a  complete  mode  structure.  (The  simple  model  is,  however,  useful 
to  gain  a  first-order  understanding  of  the  mode  structure,  as  was  done 
in  the  analysis  above.) 

An  example  of  a  complication  that  could  arise  is  with  the  "n"  mode. 
Notice  that  this  mode  can  be  expressed  either  as  E-F  or  F-E,  with  the  F 
reflection  occurring  at  either  the  transmitter  or  receiver  end  of  the 
path.  If  both  modes  propagate  and  a  tilt  exists,  the  Tq  for  each  will 
be  different;  hence,  the  mode  will  split  into  the  E-F  and  F-E  modes. 

The  degree  of  this  split  will,  of  course,  be  a  function  of  the  nature 
and  magnitude  of  the  tilt. 

Furthermore,  a  model  that  is  adequate  to  account  for  the  effects 
on  one  mode  say  a  2F  mode  will  not  be  adequate  to  explain  the  structure 
of  some  other  mode.  For  instance,  suppose  that  the  ionosphere  that 
supports  the  2F  mode  was  modeled  by  assuming  two  different  electron 
profiles  over  the  path  one  for  each  hop.  If  this  model  was  used  to 
explain  an  FEFEF  mode,  there  would  be  an  error  in  the  F  reflections, 
especially  at  the  center  of  the  path.  However,  a  model  consisting  of 
three  or  four  separate  F-layers  over  the  path  might  be  adequate  to 
explain  the  simpler  modes,  which  involve  only  two  or  at  most  three 
F-layer  reflections.  Before  approaching  such  a  modeling  technique, 
however,  the  nature  and  degree  of  tilt  over  the  path  must  be  considered. 

As  time  progresses,  the  solar  zenith  angle  changes  over  the  path. 
Correspondingly,  the  ionosphere  also  changes.  Therefore,  a  tilt  over 
the  not  i  will,  in  some  way,  be  related  to  the  time  difference  over  the 
path,  as  well  as  to  the  midpath  time.  To  calculate  the  time  difference, 
the  coordinates  of  the  transmission  end-points  must  be  known.  These 
were  found  to  be 
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Stanford:  Latitude  37.408  degrees  N 

Longitude  122.183  degrees  W 

Lubbock:  Latitude  33.600  degrees  N 

Longitude  101.840  degrees  W. 

An  approximate  value  for  the  time  difference  between  stations  may  be 
calculated  from  the  longitude  values.  Given  that  1  hr  equals  15  degrees 
of  longitude,  there  is  1.36  hr  difference  between  stations. 

The  solar  zenith-angle  variation  vs  time  of  day  for  February  was 
obtained  from  the  Central  Radio  Propagation  Laboratory  (CRPL)  Handbook 
90.  An  average  latitude  of  35  degrees  was  ass  med .  The  variation  is 
plotted  in  Fig.  52.  If  the  sun  is  not  below  the  horizon,  and  if  some¬ 
thing  is  known  about  the  nature  of  density  variation  with  solar  zenith 


Fig.  52.  SOLAR  ZENITH  ANGLE  VS  TIME  OF  DAY  FOR  FEBRUARY:  LATITUDE 
35  DEGREES. 
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angle  then  a  tilt  model  could  be  deduced.  For  example,  the  Chapman 
model, 


N  =  N  exp(l  -  z  sec  y  e  Z)  , 
o  A 

could  be  used.  This  model  could  be  initially  matched  to  the  IF  trace, 
according  to  the  midpath  time — hence,  solar  zenith  angle — conditions. 
However,  in  our  case,  the  data  used  in  the  mode  comparison  were  taken 
later  in  the  day,  when  density  conditions  are  primarily  determined  by 
electron-ion  recombination  processes.  The  sun  was  well  over  the  horizon, 
and  ion  generation  was  nearly  absent. 

The  best  idea  concerning  the  profile  changes  over  the  path  can  be 

obtained  from  the  records  themselves.  By  studying  the  properties  of  the 

IF  trace  as  time  progresses,  we  obtain  properties  of  the  ionosphere 
located  roughly  over  midpath, vs  time  of  day.  With  some  degree  of  accu- 
racy ,  one  may  then  assume  that  the  ionosphere  will  be  equivalent  at  a 
certain  time  anywhere  along  the  path.  Thus,  a  model  derived  for  the 
midpath  ionosphere  at  a  specific  time  may  be  used  later  in  the  day  to 
model  the  ionosphere  at  the  point  more  westward  along  the  path  which  has 
the  same  local  time  as  the  midpath  model. 

Now,  the  number  of  models  needed  will  depend  upon  the  degree  of 

tilt  over  the  path.  To  gain  some  idea  of  how  the  tilt  varies  with  time 

of  day,  two  very  useful  techniques  may  be  employed  with  a  sequence  of 
records  such  as  those  used  by  the  author.  For  means  of  illustration, 
these  were  tried  with  the  sequence  of  24  records  in  Figs.  1-24: 

(1)  the  IF  maximum  observed  frequency  (MOF)  was  plotted  vs  midpath 
local  time  of  day,  as  seen  in  Fig.  53; 

(2)  the  IF  virtual  height  of  reflection  was  obtained  from  the 
time-delay  information  and  plotted  vs  time  of  day  for  several 
different  frequencies.  A  few  of  these  plots  appear  in  Figs. 
54-56. 

These  two  techniques,  if  done  with  sufficiently  small  time  separation, 
and  at  enough  frequencies,  supply  all  needed  information  with  which  a 
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sequence  of  Chapman  models  may  be  derived.  However,  it  is  surely  easier 
to  obtain  the  models  by  direct  experiment-model  ionogram  comparison,  as 
was  done  earlier  in  this  report. 

At  certain  times  of  day — such  as  at  noon — the  tilts  are  very  small, 
and  therefore  do  not  warrant  the  more  complicated  modeling  analysis. 
However,  later  in  the  day,  tilts  should  be  considered,  as  seen  in  Figs. 

53  and  54-56  by  the  steeper  curve  slopes. 

We  would  need,  therefore,  a  modified  ray-tracing/interpolation 
procedure  to  accommodate  the  extra  models.  In  a  conference  with  Dr.  Croft 
the  following  ideas  were  suggested  : 


(1)  A  few  individual  models  could  be  inserted  into  the  computer. 
When  the  range  is  reached  at  which  the  ionosphere  shifts  to 
the  adjacent  model,  the  calculations  proceed  with  the  new 
model.  Figure  57  illustrates  this  technique. 


(2)  If  the  scale  height,  H  ,  and  maximum  height,  h  ,  of  the 

s  in 

Chapman  model  were  the  same  over  the  path — such  as  through  a 
sequence  of  records  resembling  Ionograms  13,  14,  and  15, --the 
change  in  the  maximum  electron  density,  N  ,  could  be  modeled 
by 


N  (6) 


N  [1  + 

°\ 


V 


k2e 


') 


where 

_  _  actual  range  -  starting  range 
1881.27  km 

and  k^,  k^  are  determined  from  the  known  N^'s  at  the  starting 
ranges  . 

This  provides  a  continuous  model,  and  necessitates  chang¬ 
ing  the  ionosphere  at  each  calculation  level . *  Similar  work 
was  done  by  Dayharsh  [Ref.  9,  p.  43]. 


*  Suggested  by  T.  A,  Croft  at  Stanford  University. 
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(3)  Higher  powers  of  6  could  be  used,  which  depends  on  the  number 
of  models  matched  along  the  path.  This  in  turn  determines  the 
number  of  starting  ranges  and  constants  k. 

If  further  investigation  of  multi-modes  is  to  be  done,  the  above 
techniques  should  surely  be  tried,  as  they  most  certainly  enable  a  closer 
match  to  nature . 

LOCAL  TIME  T,  >  T2  >  T3 


GROUND 

PATH  (km) 
TIME  (hr) 


Fig,  57.  SCHEMATIC  DIAGRAM  OF  TILT -MODELING  PROCEDURE  USING 
DISCRETE  MODELS.  Use  IF  trace  to  obtain  computer  models  at 
a  few  points  in  time  at  midpath.  Distribute  models  accord¬ 
ing  to  time  along  the  path  and  thence  synthesize  multi-modes 
from  distributed  model. 
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Finally,  as  mentioned  previously,  tilts  could  severely  hamper 

predictions  based  on  the  simple  mirror  analogy  model,  especially  over 

long  paths.  This  problem  could  be  partially  circumvented  by  examining 

records  spaced  in  time  over  the  path,  and  thence  distributing  the  F-layer 

virtual  height  over  the  path  on  a  proportional  basis.  The  h'  (virtual 

F 

height)  scale  of  Fig.  48  would  then  represent  the  F-layer  virtual  height 

at  one  end  of  the  path,  and  other  heights  would  be  related  to  this  height 

by  some  constant.  For  instance,  suppose  we  have  three  separate  layers, 

denoted  by  a,  b,  and  c.  We  might  then  have:  h'  =  1 .  lh ' ,  and  h'  =  1.2h'. 

b  a  c  a 

The  formulas  for  the  simple  layer  could  be  modified  to  include  the  new 
layers  by 

(1)  modifying  the  interpolation  ranges; 

(2)  including  more  terms  in  the  calculations  of  path  length  to 
take  into  account  the  different  path  lengths  for  each  hop. 

We  note  that  the  angle  of  reflection  at  the  ground  and  the 
top  layers,  p,  is  always  the  same. 
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IX  CONCLUSIONS 


A .  Summar\ 

The  aim  of  the  work  was  fulfilled.  An  electron  density  model  for 

the  F-layer  was  found  in  5  cases,  which  yielded  computer-synthesized 

ionogram  traces  which,  in  turn,  matched  5  experimentally  obtained  iono- 

grams  reasonably  well.  An  average  E^-layer  height  was  determined  through 

measurement  of  time  delays  associated  with  IE  and  2E  modes.  The 

s  s 

F-layer  models  and  E^-layer  heights  were  used  in  computer  ray-tracing/ 
interpolation  programs  to  produce  model  ionogram  traces  representative 
of  2F,  E-F,  FEF ,  FEFEF,  2E-F,  2F-E,  3E-F,  3F-E,  and  2F-2E  propagation 
modes.  In  almost  all  cases,  a  very  good  prediction  of  modes — based  on 
close  experiment-model  correspondence — was  possible,  as  seen  in  Figs. 
43-47.  Because  tilts  were  not  included  in  the  F-layer  models,  it  was 
not  possible  to  resolve  fine-scale  splittings  such  as  appear  in  the 
2F/FEF-E  combination. 

Successful  comparisons  to  a  simpler  mirror-analogy  model  were  made, 
which  verified  the  model  ionogram  synthesis  prediction  made  previously. 
Figure  48,  based  on  this  simpler  model,  shows  in  a  general  way  how  we 
expect  the  mode  structure  to  appear. 

Errors  were  analyzed,  and  the  degree  of  inaccuracy  caused  by  an 
improper  choice  of  Eg  height  was  derived  quantitatively  by  employing  a 
typical  F-layer  model  in  a  synthesis  procedure,  which  used  Eg  heights  of 
105,  117,  and  124  km.  An  error  on  the  order  of  ±  0.03  msec  might  be 
expected  in  an  average  model-experiment  comparison,  in  which  the  FEF  and 
E-F  modes  are  under  study. 

Tilts  were  studied,  and  it  was  suggested  that  these  could  be 
modeled  by  employing  a  sequence  of  ionograms,  to  derive  a  sequence  of 
F-layer  models,  spaced  in  local  time,  which  could  then  be  distributed 
over  the  path.  This  technique  would  enable  more  fine-scale  predictions 
to  be  made . 

B .  Applications  and  Suggestions  For  Future  Work 

Most  HF  radio  engineers  need  a  greater  understanding  of  ionospheric 
mode  structure.  A  wave  that  traverses  a  given  distance,  over  several 
different  paths,  will  arrive  at  the  receiver  with  several  different  time 
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delays,  magnitudes,  phases,  and  polarizations.  Ii,  in  addition,  there  is 
magnetoionic  splitting  into  0  and  X  rays  (which  is  usually  the  case),  one 
may  double  the  number  of  modes,  and  hence  the  complications.  A  comprehen¬ 
sive  communication  design  must  account  for  such  effects.  More  notable 
problems  would  exist  when  some  pulse  encoding  scheme  is  used,  such  as 
RTTY  or  facsimile,  where  there  could  be  severe  pulse  overlap,  which,  to 
make  problems  worse,  would  change  with  time  of  day.  Note  that  the  range 
of  modes  extends  over  some  1.5  msec;  this  number  will  be  a  function  of 
the  path  length.  Satisfactory  and  reliable  encoding  would,  under  severe 
conditions,  need  to  place  pulses  at  least  this  far  apart.  If  concurrent 
oblique  studies  are  made,  it  might  be  possible  to  make  a  judicious  choice 
of  frequencies — such  as  to  use  only  the  Eg-layer,  or  to  rely  on  the 
F-layer  near  its  MUF  where  signals  would  be  the  strongest.  D.  K.  Bailey 
[Ref.  26]  has  done  some  notable  work  in  this  area  already.  He  has 
defined  a  certain  "multipath  reduction  factor"  (MRF),  which  was  derived 
statistically  over  several  different  paths.  The  MRF  measures  the  ratio 
of  the  operating  frequency  to  the  MOF  of  the  IF  trace,  for  which  a  certain 
minimum  overlap  time  may  be  realized. 

More  recent  work  has  been  done  in  the  area  of  communication  error 
rate,  notably  by  E.  M.  Young  of  SRI  [Ref.  14],  Young  discusses  several 
sources  of  communication  error,  and  attempts  some  predictions.  His 
method  is  based  on  the  long-term  HF  propagation  predictions  of  CRPL  and 
USARPA,  and  on  estimates  of  the  "channel-scattering  function"  [Ref.  7], 
which  forms  a  basis  for  determining  time-  and  frequency-selective  fading 
effects  caused  by  ionospheric  variations. 

It  would  seem  that  (over  suitable  propagation  paths — the  shorter, 
the  better)  the  techniques  of  the  present  work  si.ould  be  of  value  in 
facilitating  communication  problem  calculations.  Through  the  use  of 
computer  ray-tracing  techniques,  with  reliable  ionospheric  models — which 
are  conveniently  described  by  exponential  functions — several  types  of 
communicative  systems  could  be  programmed.  When  a  Chapman  ionosphere 
proves  to  be  wholly  inadequate  for  modeling,  it  should  still  be  possible 
to  arrive  at  some  analytical  model  for  the  ionosphere  that  yields  well 
matched  oblique  modes. 
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Parameters  are  time  delay  and  frequency;  furthermore,  by  changing  a 
model's  parameters  with  time,  the  effective  ionospheric  doppler  shifts 
could  also  be  modeled.  Work  similar  to  this  was  done  by  Young  [Ref.  14, 
p.  30],  although  he  generally  found  poor  agreement  with  measured  data. 

Through  some  modification  of  the  computer  ray-tracing  program,  it 
is  also  possible  to  include  amplitude  variation  by  computing  range  atten¬ 
uation,  including  antenna  pattern  values,  and  estimating  ionospheric 
absorption;  Ref.  16  provides  at  least  an  introduction  to  these  computation 
techniques.  To  make  these  calculations  more  realistic,  the  earth’s 
magnetic  field  also  could  be  included  within  the  description  of  the 
ionosphere.  This  would  determine  a  "coherent  bandwidth"  through  the 
propagating  medium — as  was  discussed  by  M.  R.  Epstein  in  Ref.  27. 

Ionospheric  physicists  can  also  gain  some  understanding  of  the 
behavior  of  the  ionosphere  by  studying  the  model  variations  with  time  of 
day .  In  this  connection,  it  would  also  be  instructive  to  consult  some 
of  the  literature  regarding  the  analysis  of  ionograms  for  electron  densi¬ 
ty  profiles;  the  oblique  ionograms  could  be  transformed  to  equivalent 
vertical  ionograms,  if  tilts  were  not  too  extensive.  A  very  good  source 
is  Ref.  28,  which  is  a  special  issue  on  this  subject. 

The  types  of  modes  that  occur  are  also  functions  of  ionospheric 
parameters,  such  as  the  type  and  extent  of  sporadic-E-layers .  Further¬ 
more,  there  are  some  interesting  changes  in  structure  along  the  path 
which  occur  quite  rapidly — not  in  accord  with  normal  diurnal  changes. 

There  seem  to  be  indications  that  these  are  some  type  of  traveling 
irregularity.  If  so,  their  extent,  position,  and  velocity  will  affect 
all  the  modes  in  a  different  way.  It  seems  that  information  could  be 
gained  about  both  the  mode  structure  and  the  nature  of  the  irregularity 
if  such  a  study  were  undertaken. 

As  mentioned,  the  influence  of  tilts  should  be  explored,  and  a 
subsequent  attempt  should  be  made  at  the  fine-scale  resolution  of  nearby 
modes . 

An  exhaustive  analysis  of  multi-modes  over  several  paths  could  be 
conducted  ultimately  to  yield  some  statistical  information  regarding  the 
occurrence  of  each  type  of  mode,  and  its  relationship  to  other  types,  as 
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was  done  by  Dayharsh  [Ref.  12].  The  results,  no  doubt,  would  follow 

some  diurnal  and  seasonal  variation.  One  may  conclude  even  now  (see 

Figs.  1-24  and  the  subsequent  matches  of  Figs.  43-47)  that  earlier  in 

the  day,  the  FEF  and  E-F  modes  are  the  most  likely  candidates  of  the 

multi-mode  family.  Later  in  the  day,  as  D-layer  absorption  apparently 

drops,  other  modes  appear,  which  require  more  hops,  and  E  -layer  traces 

s 

become  much  more  pronounced.  In  general,  it  appears  that  most  probable 
modes  are  those  which  require  the  least  number  of  reflections  and  which 
are  likely  to  suffer  the  least  amount  of  attenuation. 

It  is  clear  that  reliable  prediction  of  complex  mode  structure  will 
not  come  until  E,  F,  and  particularly  Eg  structure  and  statistics  are 
understood  more  thoroughly . 

The  techniques  of  mode  estimation  (modeling)  are  so  far  a  fairly 
time-consuming  process,  and  several  more  sophisticated  computer  techniques 
would  be  needed  (as  well  as  an  on-site  computer)  to  enable  nearly  real¬ 
time  mode  analysis  and  subsequent  prediction  of  communication  effects. 
Simplified  graphical  techniques  are  still  most  efficient  for  rapid  mode 
analysis  work;  however,  there  is  still  need  for  a  reliable  graphical 
technique  with  which  one  could  quickly  explain  complex  mode  structure 
appearing  on  oblique  ionogram  data.  The  author  is  currently  working  on 
such  a  technique. 
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(Abstract  Cont) 

A  further  extension  of  the  modeling  technique 
would  include  the  insertion  of  additional  elec¬ 
trons  in  the  ionospheric  profile — E- ,  F-^-,  and  F2- 
layers — as  becomes  necessary.  Horizontal  tilts  in 
electron  density  could  also  be  modeled  by  use  of  a 
slowly-varying  profile  or  by  insertion  of  two  or 
more  different  discrete  models  along  the  propaga¬ 
tion  path. 


